Thomas-Fermi versus one- and two-dimensional regimes of a trapped dipolar 

Bose-Einstein condensate 
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We derive the criteria for the Thomas-Fermi regime of a dipolar Bose-Einstein condensate in 
cigar, pancake and spherical geometries. This also naturally gives the criteria for the mean-field 
one- and two-dimensional regimes. Our predictions, including the Thomas-Fermi density profiles, 
are shown to be in excellent agreement with numerical solutions. Importantly, the anisotropy of the 
interactions has a profound effect on the Thomas-Fermi/low-dimensional criteria. 
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Introduction. A Bose-Einstein condensate (BEG) is 
said to be in the Thomas-Fermi (TF) regime when the 
interaction energy dominates the zero-point energy [l|, \^ . 
Despite being a dilute gas, in this regime the BEG is 
strongly affected by interactions and behaves in a highly 
non-ideal manner. Furthermore, the TF description is 
relatively simple compared to full Gross-Pitaevskii the- 
ory and facilitates exact results [l]. Gonsider a BEG con- 
taining N atoms of mass m confined to a spherical har- 
monic trap V{r) = muj'^r'^/2, with harmonic oscillator 
(HO) length Zho = \/W(rnuj): and repulsive s-wave inter- 
actions with positive scattering length a. The TF regime 
occurs for large values of the parameter Na/lhoi and thus 
typically arises for dense strongly repulsive BEGs. 

Recently, a dipolar BEC was created with atomic 
dipoles polarized in a common direction by an external 
field Q. Far from any scattering resonances the inter- 
action can be represented by a pseudo-potential which 
includes the bare dipole-dipole interaction [J, Q , 
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The first term describes s-wave scattering arising primar- 
ily from van der Waals interactions, where g = Airh^a/m. 
The second term describes the long-range part of the in- 
teraction between dipoles aligned along the unit vector 
e, where Cdd is the dipolar coupling strength. An impor- 
tant quantity is the ratio between the dipolar and s-wave 
interactions e^d — Cdd/'ig [6|. For magnetic dipoles with 
moment d, Cdd = fJ'od^, where ^q is the permeability of 
free space. The signs and magnitudes of Cdd and g can 
be controlled by rotation of the polarization axis [6| , and 
by a Feshbach resonance [3], respectively. A large range 
of £dd is therefore accessible and has already begun to be 
probed experimentally 0, H, Q . 

In contrast to the van der Waals interaction, the dipo- 
lar interaction is long-range and anisotropic. This has 
dramatic consequences upon the behavior of a BEG and 
can even lead to collapse when N and/or Cdd exceed crit- 
ical values 0, [3, [ly, [ill, [I4I • Many properties of a dipo- 
lar BEG in the TF regime have already been discussed, 
e.g., its density profil e 1131 . e xpansion dynamics |13j . ex- 
citation frequencies [ij, [lJ|, rotation [l^ and vortices 
[16|. However, a systematic discussion of the criteria for 



the Thomas-Fermi regime of a dipolar BEG is currently 
lacking. In this paper we derive these criteria for the im- 
portant cases of cigar-, pancake-, and spherically-shaped 
ground states. Our approach also reveals the criteria for 
the mean- field ID and 2D regimes. 

Theory. At zero temperature the mean-field conden- 
sate wave function ^ = ^(r, t) obeys the Gross-Pitaevskii 
equation (GPE) [l|. Stationary solutions, for which ijj is 
real, satisfy the time-independent GPE, 

[-ih^/2m)V^ + y(r) + giP^ + $dd(r, t)] ^ = /it/-, (2) 

where /i is the chemical potential and the atomic density 
n(r) = '0(r)^ is normalised via J n(r)(i^r = TV. 

We assume that confinement of the gas is provided 
by a cylindrically-symmetric harmonic trap V{y) = 
^m{uj^p^ + uj'^z'^), where uix and uiz are the radial and 
axial trap frequencies, and p = \/ x^ -V y^ . The s-wave 



interactions introduce a local mean-field potential gi^P' 
while the dipolar interactions introduce a non-local po- 
tential <i>dd given by ^], 



$dd(r)=/'dV[/dd(r-r')V(r')' 



(3) 



where C/dd(r) is the second term in Eq. ([T|). 

In a time-dependent situation, where "0 is complex, ki- 
netic energy arises from both zero-point motion (density 
gradients) and velocities (phase gradients). However, in 
Eq. ([21), only zero-point energy contributes. When this is 
negligible we enter the TF regime and Eq. ([2]) becomes, 



gn{x) + m{ujip' + wf z^)/2 + $dd(r) - p. 



(4) 



Equation @ is satisfied by the well-known inverted 
parabola solution even in the presence of dipolar interac- 
tions [121 . This has the form. 



n{v) = no [1 - {p/Rxf - {z/R,f 



(5) 



where no = \bN /{SttRzR^) is the peak density. The 
density is zero beyond the BEG boundary which is spec- 
ified by the TF radii, Rx and Rz. Expressions for R^ 
and Rz are presented in [12i | . For purely s-wave interac- 
tions, stable TF solutions only exist for repulsive inter- 
actions {g > 0); under attractive interactions {g < 0), 



the zero-point kinetic energy is crucial to stabilise the 
BEC against collapse. Additionally, the aspect ratio 
K = Rx/Rz equals the trap ratio 7 — tOz/^x, while for 
dipolar interactions this is not generally true. 

It is useful to introduce a fictitious electrostatic poten- 
tial (j>{r) = J(ir'n(r)/(47r|r — r'l). This satisfies Poisson's 
equation V^0 = —n(r). For dipoles aligned in the z- 
direction, $dd can then be expressed as [12], 

(6) 



*dd(r) = -9£dd [352 
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Here the first term is anisotropic and long-range, while 
the second is short-range and contact-like. Note that the 
dipolar potential inside the inverted parabola ([3]) is p^ . 
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where /(k) = (\^1k^) j (\- k^)-Zk^zXqxA\^\ - kV(1 
K^f'l'^ lies in the range 1 > /(k) > —2 and k is deter 
mined by a transcendental equation [J, [l^ , 
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Cigar: 3D TF vs ID mean-field. Menotti and 
Stringari [i3| analysed the crossover between the 3D TF 
and ID mean- field regimes of a cigar-shaped s-wave BEC. 
We now extend their approach to include dipolar interac- 
tions. We begin by neglecting the axial trapping (cj^ = 0) 
and consider the BEC to be uniform along z with ID den- 
sity (number of atoms per unit length) rii . Then Eq. © 
reduces to the contact-like form. 



$dd(r) = ~g£ddn{p). 



(9) 



We can thus define an effective s-wave scattering length 
for the cigar ac = a(l — edd)- The dipoles in the cigar lie 
predominantly end-to-end such that for Cdd > (Cdd < 
0) the net dipolar interaction is attractive (repulsive). 

Introducing the dimensionless quantities p' = p/lx and 
''Ppip') — lxiPip)/^/ni, where Ix = y^h/mujx is the radial 
HO length, Eq. ^ becomes. 
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(10) 
The chemical potential /ii[acrii] is a function of ni: one 
can numerically tabulate this equation of state by solv- 
ing ()10|) for different values of ni. The kinetic terms in 
Eq. (|10p become negligible in comparison to the interac- 
tion terms when acUi ^ 1, in which case the TF radial 
density profile has the form. 
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with TF radius Rx — y^2pi[acni]/h!jj. Use of the normal- 
isation condition J ipp 2Trp'dp' = 1 leads to the relation, 



p,i[acni]/hujx = 2^acni 



(12) 



In the opposite regime of acUi <C 1, termed the 
ID mean- field regime [l|, the solution of Eq. PU)) ap- 
proximates the non-interacting radial HO state ipp = 
7r~^/^e~'' /^. Inserting this form into Eq. ^, we obtain 
the chemical potential perturbatively to be, 

Pi[acni]/hujx — i + 2acni. (13) 

We now wish to find the effect of finite axial trapping 
^2 7^ 0, for which ni(z) becomes inhomogeneous. We 
can easily estimate the correction to Eq. (|9l) u sing the 
TF result. Expansion of Eq. ([7]) as k ^ [l3| reveals 
that the leading correction is of order k^ and becomes 
negligible for k ^ 1 . Thus if the variation is sufficiently 
weak along z, Eq. (fTO|) still defines the local chemical 
potential along z, fii[acni{z)], and we can also employ 
the local density approximation. 



fii[acni{z)] 



2 2 

- muj^z 



/2 = M, 



(14) 



where fj, is the global chemical potential of Eq ([2]) . At the 
axial boundary Rz, this gives ^1 [flcrii = 0] -I- ^mw^i?^ = 
fi. Then, by defining the function pi — (/ii[aQ2J.(z)] — 
pi[acni = 0])/hLiJx, we can rewrite Eq. (fT4|) as (l7|. 



ili[a,niiz)]^iay2)[l-iz/Rzf 



(15) 



where a — IxRz/l^- We can determine ni(z) by using the 
equation of state /ii[acni] to invert Eq. (|15p . i.e., ni(z) = 
Pi [p,i[acni{z)]]/ac. For each choice of a we still need 
to know Rz , and this is obtained from the normalization 
condition J ni{z)dz = N which can be written as. 



K' 



\^\i-e) 



dC = A^a(l-£dd)j|, (16) 



where C, = z/Rz. On the right hand side we identify a 
natural dimensionless parameter. We enter the radial TF 
regime (and therefore the full 3D TF regime) when. 



Na{l-edd)lx/ll':>l- 



3D TF cigar 



(17) 



In this case we can use Eq. (|12l) to give the density profile, 



ni{z) = 



a(l 



Jdd) 



IxRz 
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(18) 



where Rz = [157Va(l - edd)^f/^;^]^/^- When A^a(l - 
£dd)lx/l1 ^ 1 we enter the ID mean- field regime. Then, 
Eq. (fT3|l leads to the density profile. 



ni{z) = 



1 



IxRz 

a(l-edd) V^ 



z 



(19) 



where i?z = [3Na{l—edd)ll/ll]^^'^- Importantly, at Sdd = 
1 the overall contact interactions vanish and the TF cigar 
criterion ([TT]) cannot be satisfied. When £dd = we 
retrieve the results of [rn |. 

We now illustrate the key dependence on edd- We have 
obtained numerically the ground states of the 3D dipo- 
lar GPE, following the procedure outlined in 18]. The 
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FIG. 1: Axial density 711(2) of a cigar BEG for various values 
of £dd according to the GPE (solid lines), 3D TF cigar pre- 
diction of Eq. (|18|) (dashed line) and ID mean field prediction 
of Eq. CHI) (dotted line). We employ NaL/l^ = 100, 7 = 0.1, 
LUx ~ 2-K X 50Hz and ^^Gr atoms. 



FIG. 2: 2D density n2(p) of a pancake BEG for various values 
of £dd according to the GPE (solid lines), TF pancake predic- 
tion of Eq. p6|) (dashed line) and 2D mean field prediction 
of Eq. ^ (dotted line). We employ Nall/lt = 10, 7 = 10, 
Ux ~ i-K X 5IIz and ^^Gr atoms. 



axial density profiles for a cigar BEC are compared in 
Fig. ([T]) for various values of Sdd, according to the GPE 
(solid line), TF prediction (dashed line) and ID mean- 
field prediction (dotted line). For weak dipolar interac- 
tions Edd = 0.1, giving Na{l - edd)lxlll > 1, the GPE 
profile is in excellent agreement with the 3D TF cigar 
prediction. For strong dipolar interactions £dd = 0.8, 
the ID mean-field regime becomes an excellent descrip- 
tion of the profile. For intermediate dipolar interactions 
Edd = 0.5 neither analytic form agrees well with the GPE 
results. In [13] observations of the density and expansion 
of a cigar dipolar BEC were compared with TF predic- 
tions. For their parameters, including edd ~ 0.16 and 
N mbx 10^, we obtain Na{l - edd)lx/ll « 50 > 1 thus 
satisfying the cigar TF criteria. 

Pancake: 3D TF vs 2D mean-field. We now consider 
a highly flattened pancake BEC and follow a similar 
methodology to the cigar BEC. We initially assume that 
the BEC has infinite radial extent and uniform 2D den- 
sity 712- Then, Poisson's equation V^(/) = — n(r) reduces 
to dl4) = —n{z), such that. 



$dd(r) = 2g£ddn{z), 



(20) 



This is, again, contact-like, giving an effective scattering 
length Op = (1 -|- 2edd)a- In the pancake the dipoles are 
predominantly side-by-side and so the net interaction is 
repulsive (attractive) for Cdd > (Cdd < 0). Introducing 
the dimensionless parameters z' = z/l^ and ipz{z') = 



y/l7/ri^ip{r), where Iz 
the GPE becomes. 



y^h/mujz is the axial HO length. 
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The axial TF regime exists when a-plzUi 3> 1. By 
normalising the corresponding TF density profile via 
J ipz'^z' = 1, we obtain the TF chemical potential, 

fJ,2[aplzn2]/lTUJz = {6Traplzn2f'^ /2. (22) 

In the opposite regime of aplzn2 <^ 1, we enter the 2D 
mean-field regime. Here ipz approximates the ground HO 



state ipz = TT" 



-1/4, 



^ /^ and the chemical potential is 



found perturbatively to be, 

M2[ap^2"-2]/^^2 = 1 + 2V27raplzn2. (23) 

We now consider introduce finite radial trapping LUr 7^ 
0. Expanding Eq. ^ for k ^ 00 '2^ shows that the 
leading correction to Eq. ^U\j is of the order k~^ and 
negligible when k 3> 1. Defining the function /i2 = 
{f^2[aplzn2{z)] — iJ.2[aplzn2 = Q])/hujz and employing the 
normalisation condition J n2{p)27rpdp = N wc obtain. 



2TTf3^ 



r' 
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(1 



s ds 



Najl + 2edd)ll 
It 



(24) 

where s — p/Rx and /3 — IzRx/lx- We therefore enter the 
axial TF regime, and hence the 3D TF regime, when. 



Na{l + 2£dd)ll/lt > 1- 3D TF pancake 



(25) 



For Edd = we obtain the same result as in ^21,]. The 
corresponding 2D TF density profile is. 
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where R^ = [15iVa(l + 2edd)^^//^]^/^ 

In the 2D mean-field regime of Na{l + 2edd)ll/lt < 1 
we obtain, 



n2ir) 



1 



2TTlza{l + 2£dd) V 2/ 



lzR.r 



l-ig|. (2') 



where R^ = [16iVa(l -f 2edd)lt/V2^lz]^/'^- 

The experiment [2,] featured a pancake BEC with 7 ~ 
10. For their parameters (w^ = 27r x 330, TV = 25, 000 
and Edd ~ 0.16) we find Na{l+2edd)il/lt — '^j suggesting 
that the initial BEC was not in the TF regime. 

In a very flat pancake, the net contact interactions 
become zero for edd = —0.5 and the TF regime can- 
not be supported. We demonstrate the role of £dd in 
Fig. [21 which presents the 2D density proflle of a pancake 
BEC. For edd = —0.3, the net interactions are very weak 



and the 2D mean-field prediction (dotted line) agrees well 
with the GPE profile (solid line). Conversely, for e^d = 1 
the combined interactions are relatively strong giving ex- 
cellent agreement with the 3D TF prediction (dashed 
line). For an intermediate value of edd = 0.1, the ex- 
act solution lies in between the two analytic predictions. 
Approximately spherical BEC. Away from the cigar 
and pancake limits $dd no longer reduces to a local in- 
teraction and we must consider its full form. In order to 
proceed we shall assume the validity of the 3D TF so- 
lution ([5]) and use energetic arguments to check its self- 
consistency. The combined s-wave and dipolar interac- 
tion energy for the TF solution is [13 |. 



E,= 



^^jj^ + $ddV^' 



287rm ^ 



/(K)edd|28) 



We observe that when Sdd = ^/f{i^) the net interac- 
tions vanish and the TF regime cannot be achieved. 
The kinetic energy of a spherical TF distribution is 
_Ek ~ 5Nh'^/2mRl, up to a logarithmic correction [l|. 
The TF approximation is valid when Ei/Ey^ ^ 1 which, 
using the above expressions and dropping the numerical 
prefactors (w 1), becomes {Nan/R^) [1 — /(K)£dd] ^ 1- 
From Eq. ^ we find that a spherical (k = 1) BEC oc- 
curs when Edd = £dd where e^^j = (5/2)(7^ — 1)(7^ -I- 2). 
Introducing the parameter 5 = £dd — ^dd ^^-^i expanding 
Eq. ([5]) to first order in 5 and (k — 1), we obtain the 
relation k = 1 - {75/\b)ia'^ + 2fl{\2-{^ "274-3). Fur- 
thermore, we expand the expression for R^ in [12l | and 
/(k) [22]. To first order in 5, the TF criterion for an 
approximately spherical dipolar BEC is then. 



Na\ 
ax J 



4/5 



45 



1/5 



1472(2-^7^)5 



15(1272-3-274) 



>1. 

(29) 

Since we employ R^ as the length scale rather than a^, 
this criteria has a different form from the usually-quoted 
Na/ax- We have confirmed that this gives a very good 
approximation for small deviations, typically up to 10%, 
from K = 1. It is important to note that Eq. ([5]) has 
both stable and unstable static solutions [ij]. Stable 
solutions of the full transcendental equation ([5]) for k = 1 
exist only when (127^ — 3 — 27^) > or equivalcntly 
0.5 ^7 < 2.4. Hence the factor multiplying i5 in Eq. 
is always positive and finite for the cases of interest. 



Equation p9)) reveals the sensitivity of the interactions 
to deviations from a perfectly spherical shape, as charac- 
terised by S. For a perfectly spherical BEC ((5 = 0), the 
dipolar energy is zero and the TF criterion (|29p reduces 
to its s-wave form. $dd itself does not vanish but takes 
on the saddle-shaped form <i>dd — 27inq£dd (p^ — 2z^)/5R't 
as obtained by taking the limit k ^ 1 of Eq. ^ |22| . For 
(5 > (i5 < 0), the BEC is slightly elongated (flattened) 
and the net dipolar interactions are attractive (repulsive). 
In [3|, an approximately spherical dipolar BEC was cre- 
ated {lOx « 27r X 480Hz, N = 30, 000, k « 1 and a w 5nm) 
for which the left side of Eq. (P5)) equals 45, confirming 
that it was in the TF regime. 

The TF criteria given in this paper do not specify 
whether the putative ground state actually exists. For 
a pancake dipolar BEC, radial density wave structures 
have been predicted for Sdd -^ 00 23] , and so the assump- 
tion of homogeneous radial density leading to Eq. ((20|) 
would not hold. Furthermore, a stable ground state may 
not exist; the Bogoliubov spectrum for a uniform dipolar 
BEC [10] predicts an instability to density fluctuations 
when Edd is outside of the range —1/2 < edd < 1 for 
a > 0. Although trapping can significantly extend this 
range of stability, it must be mapped out by solving the 
Bogoliubov de Gennes equations numerically [23] . How- 
ever, within this range of stability and away from density 
wave structures, the TF regime is both stable and has the 
inverted parabola profile, and so the criteria should hold. 

Conclusions. For cigar-shaped and pancake-shaped 
condensates, the non-trivial dipolar interactions reduce 
to a simple contact-like form. We identify their TF cri- 
teria, Eqns (fT7|) and ([25l) . which also determine the ID 
and 2D mean-field regimes. For the cigar (pancake) the 
net interactions are proportional to 1 — Sdd (1 + 2edd) 
such that that it becomes increasingly difficult to remain 
in the TF regime as £dd — > 1 (£dd ~^ "1/2). This high- 
lights the profound effect of the anisotropic dipolar in- 
teractions. Our predictions are in excellent agreement 
with full numerical solutions. Furthermore, for the more 
complicated case of a spherical condensate we determine 
the self-consistency criterion (|29|) for Thomas-Fermi be- 
haviour. These criteria are of relevance to current theo- 
retical and experimental studies of dipolar condensates. 
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